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Instead of using a three dimensional analysis on quasi-local horizons, we adopt a four dimensional asymp- 
totic expansion analysis to study the next order contributions from the nonlinearity of general relativity. From 
the similarity between null infinity and horizons, the proper reference frames are chosen from the compatible 
constant spinors for an observer to measure the energy-momentum and flux near quasi-local horizons. In partic- 
ular, we focus on the similarity of Bondi-Sachs gravitational radiation for the quasi-local horizons and compare 
our results to Ashtekar-Kirshnan flux formular. The quasi-local energy momentum and flux of generic isolated 
horizons and non-rotating dynamical horizons are discussed in this paper. 
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I. INTRODUCTION 

The boundary of a black hole is defined as a region of no es- 
cape of null rays, which leads to an event horizon definition. 
However, the definition of the event horizon cannot give a re- 
alistic description of how a black hole grows since it is too 
global. The event horizon can be located only after an ob- 
server knows the global structure of space-time. Hence, the 
generalization of the event horizon to a quasi-local definition 
of horizons provides an observer a possibility to detect hori- 
zons. Ashtekar et al JT] proposed a notion of isolated hori- 
zons, which is quasi-locally defined, to describe equilibrium 
states of black holes. Unlike Killing horizons or stationary 
event horizons, they do not require any space-time Killing 
fields or exclude radiations outside horizons. If gravitational 
collapse occurs, the final stage of black holes is isolated and 
in the equilibrium state therefore will not radiate any more. 
However, there should still have gravitational or matter field 
radiations outside black holes. The definitions of generic iso- 
lated horizons yield a theoretical framework to describe these 
situations. 

In order to clarify the idea of 'isolation', Ashtekar et al (Tj 
started from the most general definition of isolated horizons 
called non-expanding horizons (NEHs). It requires the in- 
trinsic (degenerate) metric q a f, on horizons to be time inde- 
pendent, i.e., £^g a f,=0, where = represents equal on horizon. 
Hence, £, which is a null normal to horizon, can be consid- 
ered as a Killing vector field of intrinsic horizon geometry. 
If one further requires the extrinsic curvature (the rotation 1- 
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form) to be time independent, it then leads to the definition 
of weakly isolated horizons (WIHs). In WIHs, the black hole 
zeroth and first laws holds. By using a freedom of rescaling 
I, it has been verified that there exists a specific I which can 
reduce NEHs to WIHs. The most restricted generic isolated 
horizons are called isolated horizons (IHs), which require the 
full derivative operator T> induced by space-time connection 
V to be time independent. 

It is expected that black holes are rarely in equilibrium in Na- 
ture. By using generic isolated horizons as a basis, the ideas 
can be generalized to dynamical horizons (DHs) definition by 
considering a space-like hypersurface rather than a null hy- 
persurface in the NEH definition. The horizon geometry of 
dynamical horizon is time dependent and it allows a quanti- 
tative relation between the growth of the horizon area and the 
flux of energy and angular momentum across it @. Ref @ 
adopted a 2+1 decomposition on a three dimensional space- 
like surface and the Cauchy data on the DHs must satisfy the 
scalar and vector constraints. Moreover, the relations between 
changes of the horizon area and energy fluxes cross the DHs 
were obtained and these fluxes has also been proved to be pos- 
itive. However, this kind of approach does not tell us what is 
the gravitational free data near horizon when considering the 
full four dimensional space-time. 

In contrast to Ashtekar et al's fl] 13 three dimensional anal- 
ysis, our works use 4-dimensional asymptotic expansions to 
study the neighborhoods of generic IHs and DHs. Since 
asymptotic expansion has been used to study gravitational ra- 
diations near the null infinity [8], it offers a useful scheme to 
analyze gravitational radiations approaching another bound- 
ary of space-time, horizons. We first set up a null frame 
with certain gauge choices near quasi-local horizons and then 
expand Newman-Penrose (NP) coefficients, Weyl, and Ricci 
curvature with respect to radius. Their fall-off can be deter- 
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mined from NP equations, Bianchi equations, and exact solu- 
tions, e.g., the Kerr solution (see Appendix [A| and the Vaidya 
solution (see Appendix [B]). This approach allows one to see 
the next order contributions from the nonlinearity of the full 
theory for the quasi-local horizons. The asymptotic expan- 
sions for the null infinity and the horizon are quite different 
geometrically. As we approach to null infinity, we consider 
an asymptotically flat space-time, however, the approach near 
horizons is not necessarily asymptotically flat. For the null 
infinity, we take the incoming null vector n a as a generator of 
null infinity to generate different cuts with respect to different 
times, say u, and the outgoing null vector I a is parameterized 
by using affine parameter r. Moreover, £ a is chosen as the gra- 
dient of the surface of a constant time u. For generic IHs, we 
take the outgoing null vector £ a as the generator of horizons 
to generate different cross sections with respect to different 
times, say v. The ingoing null vector n a is chosen as a tan- 
gent vector of ingoing null rays and also being parameterized 
by affine parameter r. The n a is the gradient of the surface of 
a constant time v. 

The main purpose of this paper is to calculate the amount 
of mass-energy fluxes crossing or near generic IHs and non- 
rotating DHs. Although there is no well-defined gravitational 
energy density in general relativity (GR), it does have well- 
defined mass or energy associated to a closed two-surface, i.e., 
quasi-locally. Without using Hamiltonian or 2+1 decomposi- 
tion on these quasi-local horizons, we can use a quasi-local 
energy formula based on spinor fields to define the mass of 
a black hole. How can one measure the quasi-local energy- 
momentum when she or he is near black holes? Our strategy 
is to use the Nester-Witten 2-form together with similar con- 
cept of Bramson's frame alignment H. The Nester-Witten 2- 
form is commonly used to prove the positivity of energy and 
can be defined quasi-locally. It has been proved to be positive 
either for the stationary space-time or the space-time with ra- 
diation, i.e, the ADM mass or Bondi mass. Therefore we use 
it for a near horizon region which can allow radiation. How- 
ever, what is a good reference frame for a strong gravitating 
field such as a black hole? To tackle this problem, we adopt 
the notion of frame alignment to find asymptotically constant 
spinors, which was done by Bramson for null infinity, and ap- 
ply this framework on horizons. After obtaining the compat- 
ible constant spinors, we use them to define the quasi-local 
energy-momentum and flux near horizons. Since we know 
the changes of the constant spinors with respect to time, it is 
possible for us to see how does the energy flux cross horizons. 

The plan of this paper is as follows. Section [II] gives a re- 
view of generic IHs and DHs. Their definitions and properties 
are discussed in terms of the NP coefficients. In Section [Til] 
we present a detail calculation of asymptotic expansions for 
IHs and non-rotating DHs. The asymptotic expansion pro- 
vides a useful scheme to analyze the space-time geometry not 
only on quasi-lcoal horizons but also their neighborhoods. In 
Section [TV] we apply the Bramson's frame alignment to IHs 
and non-rotating DHs and then find the compatible constant 
spinors. Thus, the good reference frame can be constructed 
from these constant spinors. On each cross section of quasi- 



local horizons, we find that constant spinors will satisfy the 
Dogan-Mason's holomorphic condition. In Section[V] we ob- 
tain the quasi-local energy and flux of the generic IHs and 
non-rotating DHs by using these constant spinors. Although 
our approaches, which is based on gauge fixing quasi-local 
energy-momentum and a time related condition, are com- 
pletely different from Ashtekar et a/'s, it turns out to yield 
the same result as Ashtekar-Kirshnan gravitational fluxes for 
non-rotating DHs. 

In this paper, we adopt the same notation as in the refer- 
ences JU [21 for describing generic IHs and DHs. However, 

we choose the different convention (H ), which is a 

standard convention for the NP formalism [7| (also see [5]). 
The Weyl tensor is completely specified by the five complex 
scalars ^o, ...,^4 and the ten components of the Ricci tensor 
are defined in terms of the four real and three complex scalars 
<&oo, ( J ) 22- Their definitions can be found in p. 43-p. 44 of 
|5|. The necessary equations, i.e., commutation relations, NP 
equations and Bianchi identities, for asymptotical expansion 
analysis can be found in p. 45-p. 51 of 0. 



II. DEFINITIONS OF QUASI-LOCAL HORIZONS 

In this section, we give a review of generic IHs and DHs pro- 
posed by Ashtekar et al Q]|2). The definitions of the quasi- 
local horizons yields several conditions on NP coefficients and 
these conditions are called gauge conditions throughout this 
paper. 



A. The generic isolated horizons 

We start from a 4-dimensional space-time manifold [yM ', g) 
with a 3-dimensional, null sub-manifold (A,q). The defini- 
tion of non-expanding horizon is given as follows (l): 

a. Definition. A is called a non-expanding horizon (NEH) 
if (1) A is diffeomorphic to the product SxR where S is a 
space-like 2-sphere. (2) The expansion of any null normal 
I to A vanishes, where the expansion is defined by — 
\q ah V a £b with q a b the degenerate intrinsic metric on A. (3) 
Einstein field equations hold on A and the stress-energy tensor 
T a b is such that T££ b is causal and future-directed on A. 

The geometry of A is characterized by the intrinsic met- 
ric q and the induced derivative T>, i.e., T> a — Vq, on 
A. 1 The intrinsic degenerate metric q a b on A has signa- 
ture (0,—,—). The vectors (£ a ,m a ,m a ) span the tangent 
space to A with the dual co-frame given by the pull backs 
of (n a ,m a ,rn a ). The expansions of outgoing and incoming 
null rays are defined by := ^q ab \7 a £ b = — Re[p], and 
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Q(n) == \(f b ^an b = Re [p], where <f b := -m a m b -m a m b 
on the tangent space of A. The twist and shear on A are de- 
fined as u% wist := \q a c q b d V[ 4]V [c ^ and \a shear \ := 
[\q a c q b dV (a 4) - e^]5, respectively. Since £ is the 

null normal of the null hypersurface, it implies the twist free. 
Moreover, the shear vanishes by using Raychaudhuri equation 
2 and the dominate energy condition. Therefore, the gauge 
conditions on NEH are 

K=0,cr=0,p=0. (1) 

From dTk, there must exist a natural connection one-form u> :— 
u> a dx a on A which can be obtained by 

V a i h =Lj a i b . (2) 

The surface gravity Km is defined as 

■= uj a (3) 

on NEH A (measured by £). Note that we do not have a unique 
normalization for £. Under the scale transformation £ i— > fl, 
we have u> i— ► u> + d In / and urn i— > + In / which 
leaves Eq. Q and Eq. ([3]) invariant. 

From pi, we get 

Ciq ab =Cig ab =q cb uj a £ c + q ac u b £ c = (4) 

for any null normal £ to A. In fact, I is an asymptotic Killing 
vector field as we approach the horizon even though the space- 
time metric g ab may not admit a Killing vector field in the 
neighborhood of A. 

The energy condition in the third point of NEH definition then 
further implies that R ab £ b is proportional to l a \2\, that is 
R ab £ a X b =Q, for any vector field X tangent to A. We then 
have 

*oo=*oi=*io=0. (5) 

Since £ is expansion and shear-free, it must lie along one of 
the principal null directions of the Weyl tensor. From equation 
(b) and (k) in p. 46 of (5), we have 

#o=*i=0, (6) 

so ^2 is gauge invariant, i.e., independent of the choice of the 
null-tetrad vectors (n, m, m), on A. For NEH, we also have 

dw=2(Im[* 2 ]) 2 e (7) 

where 2 e is an area 2-form. The 2-form duu can also be written 

as 

2£> [a cj b] = 2(57r - dif)m[ a m b ] . (8) 



2 For the outgoing null geodesic I, the Raychaudhuri equation yields 
= _ _ & shear & shear + ^twist + K (e)@(£) — *00- 



Im[^ , 2] plays a role of gravitational contributions to the 
angular-momentum at A. Ashtekar et al calls ui the rotational 
1-form potential and Ini[^ , 2 ] the rotational curvature scalar. 
Therefore, tt vanishes for a non-rotating A. 

Using the Cartan identity C v = di v + i v d and dTJ, the Lie 
derivative of co with respect to I gives 

£ € w a =2Im(tf 2 ) t h \a + V a (£ b io b ) =V a n {l) . (9) 

On NEH, the surface gravity may not be constant. To obtain 
the black hole zeroth law such that the surface gravity is con- 
stant, one may need a further condition, i.e., Czuu a — 0, on 
NEH. It motivates the definition of weakly isolated horizon. 

b. Definition. A weakly isolated horizon (WIH) is a NEH 
with an equivalence class of null normals under constant trans- 
formation. The flow of £ preserves the rotation 1-form ui 

C e u a =0,i.e.,[£i,V]£=0. 

From (|9j, the condition of WIH basically preserves the black 
hole zeroth law. Because £ is tangent to A, the evolution equa- 
tion is in fact a constraint. See (B21) and (B22) of [ 1 1 . There- 
fore, given a NEH, we can select a canonical [£] by requiring 
( A, [£]) to be a WIH satisfying 

C e p=0 or (i=0. (10) 

A generically admits a unique [I] such that the incoming ex- 
pansion is time independent. This result will establish that a 
generic NEH admits a unique [£] such that (A, [£}) is a WIH 
on which the incoming expansion p is time independent. 

c. Definition. A weakly isolated horizon (A, [£]) is said to 
be isolated horizon (IH) if [Cg, T>]V=0, for all vector fields V 
tangential to A and all £ G [£]. 

From this definition, we have [£g,T)]£=0 and [Ci,T>\n=Q. 
The first one gives the surface gravity to be a constant by pre- 
vious argument. So e=0. The second one gives 7r=/i=A=0. 



B. The dynamical horizon 

The generic IHs are taken as the equilibrium state of the DHs. 
The DHs can be foliated by marginally trapped surface S. 
Therefore, the expansion of the outgoing tetrad vanishes. 

d. Definition A smooth, 3 dimensional, space-like sub- 
manifold H of space-time is said to be a dynamical horizon 
(DH) if it can be foliated by a family of closed 2-manifold 
such that: (1) on each leaf, S, the expansion 6m of one null 
normal £ a vanishes, (2) the expansion 6( n ) of the other null 
normal n a is negative. 

Here, m := \ {2) q ab V a £ b and 9 (n) := \ ^q ab V a n b where 
(2)q<i b _ _( m a^j fc _|_ m b rfi a ) is intrinsic to the cross section 
S of H. From this definition, it basically tells us a dynami- 
cal horizon is a space-like hypersurface which is foliated by 
closed, marginally trapped two surface. The requirement of 
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the expansion of the incoming null normal is strictly negative 
since we want to study a black hole (future horizon) rather 
than a white hole. Also, it implies 



Re[p]=0. 



(11) 



The angular momentum associated with cross-section S is 
Jg = — <f s K ab i(j a R b dS, where ip a is the rotational vec- 
tor field and need not be an axial Killing field and K ab is 
the extrinsic three curvature J2J- ip a can be written as linear 
combinations of m a and m a , i.e., %l) a = Em a + Ern a where 
E ^ 0. The extrinsic three curvature K ab can be decomposed 
into Qab + Sab + ZW( a R b ) + BR a R b which has been ex- 
pressed in terms of Newman-Penrose coefficients |9l[l0]. If 
we choose £ to be geodesic, n a — V a u and m, m tangent to 
the two surface [9, 10 1, we get 



7^ - 
j s — 
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Re[irE]dS. 



(12) 



8tt j s 

Hence, if 7r=0, it implies a non-rotating DH. 

If we use the gauge conditions k=7t — t=tt — (a + /3)=0 on a 
rotating DH, the total flux of Ashtekar-Krishnan J2) Jl0| then 
becomes 



1 

47T 



[\a\ 2 +\7r\ 2 + ^ Q0 ]Nd 3 V 



where the gravitational flux is 



F grav = — N(\a\ 2 + \ir\ 2 )d d V. 

47r J AH 

and the matter flux of Vaidya solution is 



(13) 



(14) 



F m 



H 



T ab T a £ b Nd 3 V = — / $00^^(15) 

47T 



where we use 4nT ab £ a £ b = $ 00 . 



III. ASYMPTOTIC EXPANSIONS NEAR QUASI-LOCAL 
HORIZONS 



A. Near generic isolated horizons: vacuum 



Frame setting, gauge choice and gauge conditions 

We choose the incoming null tetrad n a — V a v to be gradient 
of the null hypersurface v = const, and it gives g ab v :a v :a = 
0. We further choose m, m tangent to the two surface. These 
gauge choices lead to 



v = /i — /i = p — p = 7 + 7 = 7r — a — /3 = 0, 



(16) 



From the definition of NEH, the gauge conditions from eq ([T]i 
can be expanded with respect to r' = r — r& as 

K = nor' + 0(r' 2 ),p = p Q r' + 0(r ,2 ),cr = a Q r' + 0(r' 2 ), 



where p := - 5 7f + 7r 7fo], a := [-So^o + ^0^0] 
from asymptotic expansions for NP equations (q) and (p) in 
p. 47 of [5 1 3 . Further we choose £ to be flag planes parallel 
on NEH, therefore it leads to 



0(0. 



(17) 



The rest of NP coefficients are 0(1). The Weyl tensor has the 
fall off (refer to equation ([6])) 



* =O(r'),*i = O(r'). 



(18) 



In order to preserve orthogonal relation £ a n a = 1, m a m a — 
— 1, £ a m a = n a m a = 0, we can choose the tetrad as 

r = (l,U,X»,X% n a = (0,-1,0,0), m Q = (0,0,£ 3 ,£ 4 ), 

where U, X 3 , X 4 are real and £ 3 , £ 4 are complex. 

We first expand NP spin coefficients, tetrad components 
U, X k , £ k and Weyl spinors with respect to r' and substi- 
tute them into NP equations which include the commutation 
relation, Ricci identities equations and the Bianchi identities 
equations. By examining equation (a) in p. 46, it implies U 
must be O(r'). X k will become O(r') since we can use the 
condition (4.86) in J9) to perform the coordinate transforma- 
tion to make X ok = 0. 

The radial NP equations 

From (n),(j,(r),(o),(q),(p),(i),(f) in p. 46-47 : 



(i = [i + (p 2 + \ \ y + 0(r' 2 ), 

A = A + (2/!oA +4Ao7o + 1'V + 0(r' 2 ), 

a = a + (X (W + /3o) + 27 a + a Q p + f§ - P V 7oK 
+0(r' 2 ), 

C 

P = A) + (Mo^o - A)( 2 7o - Mo) + a oA - P V 7oK 
+0{r' 2 ), 

p = [*°-a 7f + TO]r' + O(r' 2 ), 

a = (-P V Wo + 2l3 W )r' + 0(r' 2 ), 

7T = 7T + [2(^ + 7o )7r + 2Ao7fo + *3]r' + 0(r' 2 ), 

e = e + [2a 7To + 2/3o7i"o - 27 £o + t^o^o + ~ 7o] r ' 

+0(r' 2 ) 
K = K r' + O(r' 2 ) 

From (305) and (303) in p. 45: 

e fc = £ fe0 + O(r'), U = 2e r' + 0(r' 2 ), 



X k = 2(ir £ k0 + W t )r' + O(r' 2 ). 
From (e), (f),(g) and (h) in p. 46: 



(19) 



*o = i(-9o*i+47fo*?-3(To^)/ 2 + 0(r' 3 ), 
*i = (-5 *o + 37f *§)r' + 0(/ 2 ), 
* 2 = ^ + (-So*3 + 27fo*§ + 3/i *2)r' + 0(r' 2 ), 



*3 = *§ + (-fl *2 + 'fo*2 + 4«)*: 



0(r /: 



' From now on, we only indicate the equation and page numbers without 
mention the ref (5)- 
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The non-radial NP equations 

From (a),(b),(c), (g),(d),(e),(h),(k),(m),a) in p.46-47: 

Po = 5o K o - K 7f - k it , & =d K a , 

TTq + K = 0, 

A = 5 7r — 7r 7fo — 2A eoj 

a -- PVe O = 0, ft-PVeo = 0, 

8o/?o - Soo-o = 9oA - 9oMo = 

2Ima 7r = -2Im*§. 

(20) 

From (304) , (306) in p. 45: 

k - -2PVe , P = 0, FVlnP = ft-ao- (21) 
From (a),(b),(c), (d) in p. 49: 

*? = 2eo*? - 3/so*g, ^ = -2k *°, 



*3 



where the complex derivative is defined as V:= g§2 + « gfy, 

P(v, a; fc ) := £ 30 = -i£ 40 and P V= S . 

From (d) in p. 46 and complex conjugate of (e) in p. 46 , we 
have 

d - — c 

*o = ;r K + Po) = 2P V e = -«<>. (22) 
dv 

Using (c) in p. 46, -k = — «o- It implies k — «o = 0. 
Therefore «o is real. Using (a) in p. 46 and the fact that 
^ for the rotating NEH , we get 

k = 0, i.e., <5 £o = 0, (23) 

(i.e., P V e = 0) on rotating NEH. 
Surface gravity: from NEH to WIH 

Here we prove that the surface gravity for a rotating WIH is 



also constant. We make a coordinate choice r 



— on the 



NEH and it becomes a WIH. This gives fi = 0. Applying 
time derivative on (h) in p. 46 and using Ko = from (b) in 
p. 49, we then get eo = 0. 

It then gives us that e is constant on WIH. So the surface 
gravity re^)=Reeo is constant on WIH. For NEH, the surface 
gravity is not necessary constant. 



B. Near a non-rotating dynamical horizon: non-vacuum 
expansion 

As before, we use advanced Eddington-Finkelstein type coor- 
dinates (v, r, 9, (p). To start with, we choose n a — d a v along 
null hypersurface v = const. We have g a v. a v. a = 0. This 



gives the gauge conditions v = p—Ji = 7+7 = a+/3— If = 0. 
Further we can choose n a flag plane parallel, this implies 
7 = 0. The incoming tetrad takes the form n a = — ^ = 
g ab v.b = 62 = (0,-1,0,0). to, to are chosen to be tangent 
to the two surface. Hence, p — p, tt = T. Applying the null 
rotation type III I — > A~ x £, m 



e to, we can chose 9 to 



make Im[e] = 0. 

Gauge choices. We conclude that our total gauge choices of 
the frame in this section are 

v = p — p = p — /I = 7 = e — e = Q,ir = T = a + f3. (24) 



In order to preserve orthonormal relation, 

tn a = l,m a M a = -l,fm a = n a m a = 0, 
we can choose the tetrad as 

£° = (1,[/,X 3 ,X 4 ), n a = (0,-1,0,0), TO a = (0,0,£ 3 ,£ 4 )- 

Transfer to new coordinates (v,r',6,(f>) 

Since the dynamical horizon is not a null hypersurface, we 
have to transfer to a 'good' coordinates which can capture this 
picture. We defined a new radius coordinate by using 

r' = r — r^(v). 

To tackle the problem of the power series expansion near a 
non-null hypersurface, we make a coordinate transformation 
by using dr' = dr — r\dv to transfer coordinates into the 
new coordinates (v,r' ,6,(j>). The tetrad in the new coordi- 
nates (v, r', 8, (f>) are: 

^ = (l,(7-f A ,X 3 ,X 4 ), 

n a = (0,-1,0,0), TO a = (0,0,£ 3 ,£ 4 ). 

The metric has the form 

g ab = £ a n b + n a e b _ m a-b _ -a^ 

where it's components are 



9 



01 



1, 3 UU - 9 Uk = 0, g 11 =_-2U + 2f A , 

g lk = _ X k^ g mn = + £ m £«). 



H is a non-null hypersurface, since g ab r' a r\ b = —2U + 2f A 
in this new coordinates. 

Gauge conditions of non-rotating dynamical horizons 

On the non-rotating dynamical horizons, the null normal £ a 
is hypersurface orthogonal, therefore it is geodesic and twist- 
free. Moreover, it is the null normal of the marginal trapping 
surface, therefore it is expansion-free on the dynamical hori- 
zon. However, unlike the isolated horizon, the shear is non- 
vanishing on the dynamical horizon. Here we only consider a 
non-rotating dynamical horizon, so n vanishes on the bound- 
ary. In terms of NP, these conditions are 

K=p=ir=Q 
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on the non-rotating DHs. Therefore, the fall off of NP coeffi- 
cients and the components of the null tetrad are 

U,Z k ,X k ,ii,e,a,P,\,*,= 0(l), 
k, p, it = O(r'). 

The falloff of the Weyl tensor and stress energy tensor 

Because we want to compare with the Vaidya solution, we 
assume that the falloff of the Weyl tensors and Ricci tensors 
is the same as that of Vaidya solution. Therefore, we take the 
Weyl tensor as 

$ = O(r'), *i = O(r'), # 2 = $ 3 = *4 = 0(1) (25) 
and the Ricci tensors as 

$oo = 0(1), $22 = *n = $02 = *oi = *2i = 0(r').(26) 

We note that the coefficient of order O(l) of the tetrad com- 
ponents U, X k , i.e., Uq, X ok can be made to be vanished by 
using coordinate transformation. 

The radial equations 

From (n),(j,(r),(o),(q),(p),(i),(f),(c) in p. 46-47 : 

p = /io + (Mo 2 + A Ao)r' + 0(r' 2 ), 

A = A + (2 Ai0 A + ^)r' + O(r' 2 ), 

a = ao + (a p _-A ao + *3)f' + 0(r' 2 ), 

/3 = -a + (a \ - p a a )r' + 0(r' 2 ), 

p = (a A + ^)r' + O(r' 2 ), (27) 

a = <t q + p Q <r a r' + <3(r' 2 ), 

7T = *§r' + ( Mo ^ + Ao^o)r' 2 + 0(r' 3 ), 

e = e + $Er' + O(r' 2 ), 

K = r A ^>lr' + 0(r' 2 ), 

From (305) and (303): 

U = 2e r' + ^y 2 + 0(r' 3 ), (28) 

X k = (^fe0 + ^fe0) r /2 + O ( r /3) 

From (e), (f),(g), (h), (i),(j) and (k) in p. 49-51 : 

$o = (-3a $ 2 l + Ao$ l )r' + 0(r' 2 ), 

*i = (-2a *g-5 *»)r' + O(r' 2 ), 

$ 2 = ^ + (-5o^ + 3/z *°-a *0)r' + 0(r' 2 ), 

*3 = $g + (-9 ^+4Mo*§)r' + 0(r' 2 ), 

$4 = + + 0(r' 2 ), 

$oo = $ l + 2 Mo $oV-' + 0(r' 2 ), 

$n = ir A $V 2 + 0(r' 3 ), 

$oi = -f A $? 2 r' 2 + 0(r' 3 ), 

$22 = $^ 2 r' 2 + 0(/ 3 ),$ 12 = $? 2 r' 2 + 0(/ 3 ), 

$02 = $g 2 r' 2 + 0(/ 2 ). 

(29) 

The non-radial NP equations 

From (304) and (306) in p. 45: 

P V r A = 0, P = r A (A P + p P), a = §P V In P. 

(30) 



From (a),(b),(c), (g),(d),(e),(m),(l),(n) in p. 46-47 : 

r A ($ 2 + cr A ) = -$g - cr CT , 
cr = 2e cr + r A p a , 

k q = -2P V eo = ?"a$3> 
A = -2e A + /io^o + rA(2^ Ao + $^), 
"o = r_ A (p a - A a + |$") - ao^o, 
*3 = 9oMo - 9oA , 

$° = ^(2 V InP V InP- V InP- V InP) - A cr , 
Im$° = -Im(A CT ), 

$° = <9„/i + 2e Mo - A cto - ^a(m 2 + AoA ). 
From (a),(b),(c), (d) in p. 49: 

rA(2ao$§+ao$^)--PV$g 0l 

d v ^1 = P V $[j + 2e $? - 3«o#§ + 7fo*8o 

-2f A ($g 1 -*i) ! (31) 
5,$° = p <S>°oo + ^a(-9o$§ + 3p *° 2 - o- *2), 
d v *° 3 = 9 $« - 2e $ 3 l + r A (-9o$? + 4 Mo $§), 
0„*2 = 5 $g - 3A #g - 4e $ 4 l + r A $I- 

For example, &q has a next order contribution from asymptotic 
expansion since &q — 2eoO"o + r A <7i, where <j\ = p<tq. 



IV. CONSTANT SPINORS FOR QUASI-LOCAL 
HORIZONS 

A. Constant spinor s for the generic isolated horizons: Frame 
alignment 

In this section, we adopt a similar idea of Bramson's asymp- 
totic frame alignment for null infinity [4| and apply it to set 
up spinor frames on the quasi-local horizons. We define the 
spinor frames 

Z A ±={\ A ,n A ) (32) 

where A^ = \\oa — Ao^a, Pa — P\oa — Po^A- We expand 
Ai,A as 

Ai = A?(«, 6, 0) + X\(v, 9, <f>)r' + 0(r' 2 ), (33) 
A = Ag(«, 0, 0) + Xl(v, 9, <f>y + 0(r' 2 ). (34) 

Here Ai is type (—1, 0) and Ao is type (1,0). 

Firstly, we demand the conditions on spin frames to be paral- 
lelly transported along null normal £ a direction on the quasi- 
local horizons, so 

lim r ^ DZ A A = Q, (35) 

and also the conditions of the frames on different generators 
on the quasi-local horizons are: 

linv^o 6Z A ± = lim r ^ SZ A ± = 0. (36) 
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We then get six conditions 



We define the spin frame Za — = [Xa, /m) where 





e A 


= 0,i.e., p o A° = 0, 


(37) 


A?4 


eoA? 


= 7r Ag,i.e., p A" = 7t Aq, 


(38) 




a Aj] 


= o, 


(39) 


o A°- 


MoAq 


= o, 


(40) 






= o, 


(41) 


9oA? + 


°0 A 


= 0. 


(42) 



for the constant spinors Xa on the generic isolated horizons. 
In (42 1 and also the following, we use the symbol — a' to rep- 
resent the one of the NP coefficients, Ao (the leading order of 
NP shear of n), in order to avoid confusion with the spinor A^ 

We use the condition ( |37j > and the fact that p 5o = QoPo on 
the horizon. Apply p on ( 39 1, we find 



= p 5 A|j = 5 p A° = 0. 



(43) 



So condition ( 37 1 and ( 39 1 are compatible 



Apply b on fi0] > and use condition p7\ and ( |38| ), we have 

= <^A°. (44) 

Hence the condition ( [37| , p8j ) and pOj l are not compatible 
unless = 0. 

From the previous analysis, we conclude that the compatible 
frame alignment conditions for the generic isolated horizon 
are (P37}, d39} and (HOI. Equation d39} and (HDl are Dougan- 



Mason's holomorphic conditions |6|. Here we see that the 
conditions of the spinor field to be asymptotically constant 
on NEH implies the Dougan-Mason holomorphic conditions 
on the cuts of the NEH. These equations will be used to- 
gether with the Nester-Witten two form to define the quasi- 
local energy-momentum. The time related condition p7] i will 
tell us how the energy momentum changes with time along 
NEH and will be useful to calculate the energy flux across the 
horizon. 



B. Constant spinors for a non-rotating dynamical horizon 



\\Oa — Xqla, 



(45) 
(46) 



and 



A = X (v,e,cl>) + X 1 (v,e, ( t>y + 0(r' 2 ), (47) 
Ai = \°i(v,e,® + \l(v,9,ct>y + 0(r' 2 ), (48) 

where Ai is type (— 1, 0) and Ao is type (1,0). 

Firstly, we demand that the conditions on spin frames to be 
parallelly transported along the non-rotating dynamical hori- 
zon generators R a direction on H, so 

linv^o RZ A - = liov-oOD - r A A)Z A A =0, (49) 

where we use R a — (£ a — r&n a ). The conditions of the 
frames on different generators on H are: 



lim r /^o SZa — = lirrif'— >o 5Za — = 0. 



(50) 



Therefore, from o A (D — rAA)A J 4 in equation (49 1 we get 

(L>-r A A-e + f A 7)Ao + (K-r A T)Ai = (51) 

where we use 7 = 0, tt — r=0, n — 0(r'). The total six 
equations that include two time-related conditions and four 
spatial-related conditions are 



Aq - eoAg 

A? + eoA? 
a A° + a A? 
^oA? — ^qAq 



is' x° 
°o A o 



5 A5 — AqAq 



0,i.e., pAg = r A p'A{j = 0, (a) 

0,i.e., pA? = f A p'A? = 0, (b) 

from d\ + aXi = 0, (c) 

from 5Ai — /zAo = 0, (d) 

from 5Ao + pXi = 0, (e) 

from 5Ai + AA = 0, (/) 



From the commutation relations, we have p 5 — 5op =£r 9o- 
From (a), we have doPoAo = 0. From (c) and use the fact 
p (7o = 0, we have p 5oAg = — ooPoA?- To make this rela- 
tion compatible with commutation relation, we need an extra 
condition 



The dynamical horizon is a space-like or time-like hypersur- 
face if the horizon is expanding f A > or contracting. From 
this we know that it's generator is non-null and it is a time- 
like or space-like generator respectively. To understand this, 
we have to work in the (u, r' , 0, 0) coordinate. The dynam- 
ical horizon generator R is tangent to the dynamical horizon 
where Crv = 1 so R = 4-. We have the relation 

R a =£ a -r A n a = -^-, 
ov 

R a = 4 - r A n a = -2r A dv - dr', R a R a = -2f A . 

So if f A > then the H is a space-like hypersurface. Since 
g ab r[ a r[ b = 2f A . 



p A° = -a AEJ. (g) 

To examine whether this extra condition is compatible with 
Dougan-Mason holomorphic condition, we apply p on (d). 
From (d), we have 

p 5 A? - (rroAo + y° 2 )X° . (52) 

From (g) and together use the fact 5 ero = and commutation 
relation of [5, 5'], we have 

5 oPo A° - a 5 A? + (aoAo + (53) 

Therefore, this extra condition is compatible with Dougan- 
Mason holomorphic condition. 
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From the previous analysis, we concludes that the compati- 
ble frame alignment conditions for the non-rotating dynam- 



ical horizon would be one time related condition (54 > (i.e.. 
(a)), two Dougan-Mason holomorphic conditions (55 1 (i.e.. 



(c)) and (56 1 (i.e., (d)) and one extra condition (57 1 (i.e., (g)): 



Mo 
5 A° + a A° 

So A? — MoAq 



rAp'A° = 0, i.e.. 
0, 

°'_ 

^SoAq, 



e X° = 0, (54) 
(55) 
(56) 
(57) 



where first three conditions are similar with the compatible 
constant spinor conditions for the generic isolated horizon. If 
the shear term vanishes, we do not need this extra condition. 



Remark. The Dougan-Mason holomorphic conditions ( 55 i 



and ( 56 1 will tell us how to gauge fix the horizon quasi-local 



energy expression and the time related condition (54 1 will tell 
us how the energy momentum change with time along the dy- 
namical horizon. 



V. ENERGY-MOMENTUM AND FLUX 

A. The quasi-local energy-momentum of an isolated horizon 

By using Nester-Witten two form together with the compat- 
ible constant spinor conditions which are Dougan-Mason's 
holomorphic conditions ( [39] > and ( |40| > for the NEH, the quasi- 
local momentum integral near a NEH is 

I(r') = -i/ Sr [Ao'aAi-Ai5Ao' + Aoa'Ai/-Ai/5'Ao 
-A A '(/x-(-7l) - XiXy(p + p)]dS 



MoAqAq/ + 0(r')]dS. 



(58) 



Moreover, the horizon momentum PaA' can be written as 



Paa'(Sa) = I(r A )X A XA> 



(59) 



where A ,4 is constant spinor on two surface of NEH. From 
the result of the asymptotic expansion for the generic iso- 
lated horizons, we can re-interpret the quasi-local energy- 
momentum integral of the generic isolated horizons (NEH) as 



I(r A ) 



1 

47T 



2e 



A^A ,d5 A (60) 



where ^° = M + iL and ^o^o = A — iL with M, L, A are 
function of (v,9,(j)). We compare with Kerr solution in the 
Appendix [A] 



B. News function of the generic isolated horizons 



In order to match the Kerr solution that its flux vanishes, we 
rescale the spinor field. Firstly, the constant spinors A{j and A] 1 



are rescaled by using the following relation 

Ag = Age"/ e °^, A? = A?e" S e ° dv , (61) 
and it yields the new rescaled momentum integral 

p -2fe dv. 



I{r A ) 



' J I(r A ) = -i- I (Xo\ \ldS A .(62) 



The three compatible conditions (37i, (39i and (40 » then be- 
come 



X° =0, 6 X° = 0, 5 A° - p X° = 



(63) 



where we use <5o£o = from asymptotic expansion and they 
are still compatible under rescaling. 

By using this new rescaling constant spinor frame, we apply 
the time derivative on the quasi-local energy-momentum of 
NEH (|58]> and thus we get 



1 

47T 



poX X tdS A 



(64) 



We said that ( |64] l is quasi-local energy flux near NEH. Here 
p is related with the mass loss or gain, hence it is the news 
function of NEH. However, since the black hole laws do not 
hold on NEH, we conclude that the news function of NEH is 
not physically reasonable. This is due to that one cannot mea- 
sure the correct temperature of a black hole since he or she 
uses a bad thermometer, which corresponds to the normaliza- 
tion of I. If one choose the canonical [£] which yield WIH, 
then po vanishes. 



C. Conserved quantities of the generic isolated horizons 

We now consider the absolute conservation law that G m — 
on the generic isolated horizons. From (b) in p. 49, we have 
^2 = 0' therefore, we can find ten conserve quantities which 
are 



2^2, m 



^° 2 dS. 



(m = -2,-l....,2) 



(65) 



Here these conserved quantities corresponds to three different 
type of generic IHs. These conserved quantities are related 
with mass and angular momentum of the generic IHs. 



D. Quasi-local energy momentum of a non-rotating dynamical 
horizon 



By using the Nester-Witten two form and Dougan-Mason's 
holomorphic conditions ( f55] > and ( |56*] l which we found in per- 
vious section, the quasi-local energy-momentum integral on 
horizon can be expressed as 



I(r A ) = - 



1 

4tt 



Po^oXaidS A . 



(66) 
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Then use (h) in p. 46, the quasi-local energy-momentum in- 
tegral for a non-rotating dynamical horizon can be expressed 
as 



The energy-momentum integral (68 1 then becomes 



J(r A ) = - 



1 



1 



4tt ./ 2e, 



+rA(Vo + c'oo-'o)]AoA ,dS , A. 



(67) 



(68) 



E. Flux expression of a non-rotating dynamical horizon 



Analysis and properties of NP equations 

<1> From (a) in p. 46, because of dominate energy condition, 
we have r A ( 1 i>2 — coer'o) = — ^oo — "o^o < and 
f A > 0. It implies ^ - ow'o < 0. 

<2> From imaginary of (1) in p. 46, it implies im^ = 
((Tocr'o — oW'o)- If Im^ — f° r shear non-zero, 
then it implies cr'o = — A&q where A is real and time 
independent. 

<3> Use <2> and from the result of asymptotic expansion 
(b) in p. 46 and (g) in p. 46, we get = (-4Ae + 

Af A no + Ho)a + r A * 4 . (1) If r A = 0, it implies 
(— 4Ae + Mo) ^ = 0. So either a — or [i — 4Aeo- 
For the later case r A = does not imply shear van- 
ishes. Since r A should imply shear vanishing, this case 
does not satisfy the boundary conditions of NEHs. (2) 
If <7o = 0, then there are two situation, (a) r A = then 
it goes back to NEH. (b) *° = but f A not vanish. It 
is similar to Vaidya. 

<4> From (g) in p. 46, if cr'o = 0, i.e., A = 0, then ^o^o = 
—rA^I- (1) If f A = since /io cannot be zero then 
(Jo = 0. (2) If Co = 0, then there are two situation. (2') 
rA — then it goes back to NEH. (2") = but f a 
not vanish. It is similar to Vaidya. 

From <3> and <4>, it would be more reasonable for us 
to chose cr'o = 0, i.e., A = 0. 



—, then u is not 



<5> If we choose the coordinate r A - 

related with 9 or <f>. So 5o/io = 0. Together with the 
choice cr'o = 0, we then have ^§ = from (m) in p. 
47. 



From the above analysis and reasons, we choose r a = — ^ 
and set 

v'o = 

to derive flux expression. From (g) in p. 46, (b) in p. 49 and 
(h) in p. 46, we have 

MoCTo = -rA^t (69) 
*2 = Mo^oo + SMoeorA + Mocroffo, (70) 



I{ta) 



1 

47T 



HoXoX Q ,dS A (i) 



lf§/o^>S A (ii) 



Now we apply time derivative on the energy-momentum inte- 
gral and also we use the time related conditions ( |54| ) of spinor 
field. 



Apply time derivative on (i), we get 

I(r A ) = — j> iio\° >? ,dS A . 



Apply time derivative on (ii), we get 



1 

47T 



^[w + $l ]\p? a ,dS A . 



(72) 



(73) 



Finally, we get the quasi-local energy-momentum flux for- 
mula which is obtained from asymptotic expansion for 
the non-rotating non-spherical symmetric dynamical horizon 
(shear non-vanishing on the boundary) is 



I(r A ) = 



i I fi XoX" ,dS A 



1 



4 _ j ^[aoao + ^oCdSA (74) 

which is quite similar with Ashtekar's flux equation for dy- 
namical horizon. 

Therefore, /io is the news function of non-rotating DHs that 
determines gravitational radiation and matter field radiation. 
The shear square term | cxo 1 2 is believed to be related with 
gravitational radiation, since the gravitational radiation comes 
from the non-spherical symmetric gravitational collapse. In 
the later section, we will see more detail about how this equa- 
tion can be compared with Ashtekar's flux formula. 



F. The relationship between Ashtekar-Krishnan flux and our 
flux formula 



From (i) 

From the choice of /io = — we have = ^f- = ^ ( 2 >R 

where the two scalar curvature is = \ (The metric of a 

two sphere with radius r A is dl 2 — —r\(d0 2 + sin 8d<fi 2 ).). 
We use these relations and substitute them into the the time 
derivative of (i) (recall ( |72] i) in previous section, we then ob- 
tain 



I(r A ) 



AoAqAq/ dS A 

-RAoAq/oJS'a- 



2eoMo- 



(71) 



dJ J_ 

dv 4-7T 

- 1 I dr * (2),,^" 



8ir J dv 
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Integrate the above equation with respect to v, we then have 

dJ=i / WR\ \° ,dS A dr A . 



8tt 



(75) 



We recall that the Ashtekar's total flux formula J2) is 



From 1. and 2. and with the choice of N — 2AqAq, > 0, 
we then have — „ M ° = 2. Our formula is then the same 
as Ashtekar-Krishnan's formula (|80|. Therefore, the surface 
gravity is = . The differential of the horizon 
area is dA = 8irr A dr A and K^dA — Airdr A . For the time 
evolute vector t a = N£ a , the difference of the horizon energy 
E 1 can be expressed as 



" 'matter ~t~ Fgrav — 



4r~ f {2) RNd 3 V (76) f /" J1V , 1 /" 

i6tt y Aff J = J (7 ' a) = 4^ y 



which F matter + F grav is equal to flux dl and d 3 F = dr A dS 

on horizon. Therefore, if N = 2AqAq/, the our flux for- 
mula from equation (i) is completely the same with Ashtekar- 
Krishnan's formula |76]i. 

From (ii) 

Recall the Ashtekar-Krishnan gravitational flux expression 
([RJis 



F, 



, -^-/ iV(|a| 2 + K| 2 )d 3 F 



(77) 



where it satisfies the gauge condition of our asymptotic ex- 
pansion. The matter flux expression of Vaidya solution is 



F m 



T ab T a e b Nd 3 V 



H 
1 

4tt 



$oo Nd 3 V 



(78) 
(79) 



where we use AmT ai) i a l h = 4>g . The total Ashtekar-Krishnan 
flux of non-rotating DH [2| [ 1 1 then becomes 



iW= ^ J[\a\ 2 + <f> Q0 ]Nd 3 V. 



(80) 



Recall our formula f73] l in previous section and then integrate 
it with respect to v, we have 



dI(r A ) = -—J ^[a o a + $ u 00 ]\»T Q ,dS A dv 



47T 7 2e rA 



where d« = (Since dr' vanishes on horizon.). 



1. Using (h) in p. 46 and (1) in p. 46, we have 

1 1 

'a 



2eoMo = ~P V a - P V a Q + 4a a = = -^Mo 



where we use the fact that for a sphere metric ao 
- V l^ e . It then implies 



Mo 



-4e . 



(81) 



2. Because u is arbitrary, one can always rescale v so that 
we can chose r A = 1 . 



dI(r A ) 
7;dr A = 



1 

47T 

8tt 



[IH 2 



$oo]^ 3 ^ (82) 



(83) 



which it is always positive. We can get a generalized black 
hole first law for non-rotating DHs 



(i) 
8tt 



(84) 



VI. CONCLUSIONS 

We know that the news function determines the gravitational 
radiation near the null infinity. It is the time derivative of the 
shear of outgoing null tetrad (<f ). Similarly, the horizon news 
function can determine the gravitational radiation and energy 
loss or gain near the horizon. Our work of the asymptotic ex- 
pansions near quasi-local horizons does not require the space- 
time is asymptotically flat. Using the method of asymptotic 
expansions we find that the corresponding news functions near 
quasi-local horizons are the time derivative of the expansion of 
the incoming null normal (/to)- If we a Pply tne ti me derivative 
on the Bondi mass, it will give us a minus quadratic of time 
derivative of the shear term. It is always negativel Therefore, 
for an isolated gravitating system, it will always lost mass and 
the gravitational wave will carry gravitational radiation out to 
infinity. 

The news function of dynamical horizon determines the quasi- 
local energy flux cross the dynamical horizon to be always 
positive, if the dominate energy condition holds. Then the 
dynamical horizon will gain mass and the horizon will grow. 
In this paper, we have investigated the quasi-local energy of 
quasi-local horizons for which the physical energy conditions 
hold. 

Asymptotically constant spinors can be used to define the 
quasi-local energy-momentum of the horizon. Searching for 
the compatible conditions of constant spinors of horizons of- 
fers us a way to chose for the proper reference frames when 
measuring these quasi-local quantities. In practical, Dougan- 
Mason's holomorphic conditions refer to how to fix the gauge 
(good measurement) in quasi-local energy-momentum ex- 
pression (Nester-Witten two form) on each cross section of 
horizon. The time related condition can tell us how the quasi- 
local energy-momentum change with time along horizon. 

We find that the news function exists for the NEHs. It indi- 
cates two possibilities: one is that there are gravitational radia- 
tions outside the NEHs, and the other is that gravitational radi- 
ations cross to the equilibrium black hole but the horizon area 
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will not increase. For the first one, we cannot prove whether 
news function corresponding to gravitational radiations out- 
side the NEHs. For the second one, it would be against black 
hole thermodynamic laws. However, since the news function 
of the NEHs can be made to be vanished while we make a spe- 
cial choice of affine parameter ta = we conclude that 
this is due to a bad measurement (bad gauge choice) of gravi- 
tational flux. This result refers to that a generic NEH admits a 
unique [£] such that (A, [£]) is a WIH on which the incoming 
expansion /io is time independent HI . 

Laws of black hole dynamics 

Zeroth law For NEHs or DHs, the zeroth law will not hold. 

If we make a coordinate choice r a = — 1/Mo for NEH, then 
WIH will preserve the zeroth law mechanics. 

First law From our construction, we can derive the general- 



ized first law eq. ( 84 1 for non-rotating DHs. 



Second law It is obvious that our flux formula for a dynam- 
ical horizon is related with the black hole area law. If the 
news function [iq is positive, i.e., the positive energy flux. 
From the fact of the area of cross section changing with time 
4-dS& = ^Tj\fiodSj\, we know that the black hole area 
is always increasing if the dominate energy condition holds. 
Therefore, the black hole second law can extend to a dynami- 
cal horizon. 

For a spherical symmetric dynamical horizon, the contribu- 
tion of energy radiation is purely from the matter field. There- 
fore the gravitational contribution is from the non-spherical 
symmetric term: the shear term. Our result about the gravita- 
tional radiation of non-rotating DHs that comes from the next 
order contribution of asymptotic expansion is proportional to 
the shear square and always positive. This result agrees with 
perturbation method of energy flux cross event horizon [Done 
by Hawking and Hartle, see Chandrasekhar 0] and Ashtekar- 
Krishnan's flux formula for non-rotating DH. 

One can refer to Chapter 6 of [9] or ref [10] for the relations 
between energy fluxes and area balance of rotating dynami- 
cal horizons. The future work on gravitational radiation and 
angular momentum flux of rotating DHs is in progress ifTOl . 
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APPENDIX A: QUASI-LOCAL MASS FOR THE KERR 
SOLUTION 



tion for Kerr horizon mass Q. He used the Dougan-Mason 
mass definition and solve the holomorphic conditions to get 
the exact expressions for the constant spinors. To calculate 
the Dougan-Mason mass on the bounded two sphere, i.e., at 
the quasi-local level, we need a tetrad to satisfy the gauge con- 
ditions (p, fi are real). Therefore, we need to make a null rota- 
tion of Kerr solution in advance Eddington-Frankelstein coor- 
dinate. By using the constant spinors, the new NP coefficients 
and the Dougan-Mason mass, we can calculate the quasi-local 
mass for Kerr. In order to compare with our work, we expand 
it with repect to angular parameter a when a is small. 

To make our tetrad to be compatible with the method of the 
asymptotic expansion for Kerr, we transfer it to the Bondi co- 
ordinate. Thus we have the desire gauge conditions in Bondi 
coordinate when a is small in Kerr metric. By using the re- 
sults of the Kerr solution in the approximate Bondi frame in 
Subsection |A 1| we can solve the holomorphic conditions and 
calculate the quasi-local mass by using Dougan-Mason's def- 
inition in Subsection IA 31 



In Subsection |A4| we check the Kerr solution by using our 



quasi-local mass that obtain from asymptotic expansion ( 60 1 
and see whether it can compare with the other two resuls in 
the previous sections. Our conclusion is that all these three 
calculations are the same up to the second order when power 
expanding with respect to the angular momentum parameter 



1. Kerr solution in the approximate Bondi frame 

We transfer Kerr solution from advanced Eddinton- 
Frankelstein cooordinates to Bondi coordinate near horizon 
(detail see Chapter 6 of |9|). The Kerr metric in Bondi 
coordinate (v, r', 8, <j>) near the horizon is therefore given by 

5 00 = 0(r'% ff 03 = O(r' 2 ), g 12 = 0(r'% . 9 02 = 0(r' 2 ) 



r A 2 + a 2 r 2r A a 2 sin 2 6 a? sin 2 6»(r A - M) . . 

r 



V 2 
A 



(r A + a 2 )£ A 



,13 



+ o(/ 2 ), 

a a(r A - M) (2r A 2 + 2a 2 -a 2 sin 2 i 
~ ~£a (r A 2 + a 2 ) 2 X A 



9 22 = --^+2-^V + 0(r' 2 ), 



33 = 

y - -2 



2(r A 2 + a 2 ) 2 SA 



Sa ' "S A 2 

1 ( a 2 (2r A 2 + 2a 2 -a 2 sin 2 ( 
sin 2 6>£ A 
2a 3 sin 9 cos 6 

> A 2 + a 2 ) 2 X A 
■ rA - M V + 0(r' 2 ), 



O(r'), 



-r' + ()(r ). 

7 " = -2 r " 



Here we present three different calculations for Kerr horizon where ta := r + = M + \J M 2 — a 2 and Sa 



mass. In Subsection A 2 we re-calculate Bergqvist's calcula- 



a 2 cos 2 tt. 
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Now we consider the case of slow rotation so that a is small, 
the tetrad components in the Bondi coordinate (v, r', 8, tj>) are: 



i a = (l,[//,0,— + Dr'), 

T A 

n a = (0,-1,0,0), 



(Al) 
(A2) 



m a = -J— (0,0,1-— ,^(1-— )), (A3) 
V2r) A va sva.6 t] A 

where n A '■= r A + ia cos 9, U :— rA 7 M and D := 
a(2i A -M) ^ j^jp coe ffj C j ents an( j Weyl tensors are: 



K = u = A = z^==0, 

U(-r A + r')r' 



P 



=0, 



(VA - r')(r) A - r') 

—r A + r' ^ r A 
(VA-r')(r] A -r / ) S A ' 
_ iy/2Dr) A sin 8 „ i\f2Dr\ A sin 9 



4(t?a - r 1 ) 
iD sin 8rj A 
y/2iacos(28) 



(3 = — —iD sin 6rj A + — — 

8 42jA 



7 



4sin6>£ A ' 

U[(r' -r A ) 2 + a 2 cos 2 8 + ia cos 8r'} ; [7 

2[(r' - r A ) 2 + a 2 cos 2 0] = I" 

m cos _ 

-,7 + 7=0, 



(A4) 



2S A 
7T = a + (3, 
*o = 0,1-1=0, 

T T „ i-D COS ? 2 o 2 • 2 /,x 

Im*2 = (r A + a 2 cos 2 6> - a 2 sm 2 8) , 



#3 = 



i\/2sin 8r A ri A . ~ , 9 „. 

— , 3 ' [D£ A + 2ia 2 cos 2 0] , 



#4 = 0. 



From this Kerr tetrad in the approximate Bondi frame, the NP 
coefficients satisfy 

iy^fi — JZ==it — a — /3==7 + 7=e — e=0, 

7T=T, p=p, /i < 0. (A5) 



2. Bergqvist's calculation 

After making the null rotation, the complex spatial null tetrad 
is 



1 



rn 



-[do 



sin%| + a 2 ) 

and the NP /3 then becomes [ 3 1 

„,„, 1 . „ 2iasin8 ia(r — M) sin 9, 

/3(9) = _ cot 6> ^ — 

2\/2ry A Va (»i + a 2 ) 



The Dougan-Mason's holomorphic conditions can be written 
as two equations in component. The one is 



= 5A = (S - f3)\ Q , 



(A7) 



where the shear vanishes. The solution has the form Ao = 

C*_iA_i(6»)e-^/ 2 + CiA 1 (8)e l ' t >/ 2 . The A- X ,A X satisfy the 
differential equations 



d 8 A_ x -(G + H)A_ x =0, 
d e A 1 + (G - H)A 1 = 0. 



(A8) 
(A9) 



With the aid of using Maple, we can find the solutions are 

A (ff\ — r f G + Hd9 — sin ( 9 / 2 ) PYr ,f "tcos9(r-M-ia) >, 
— 1 \ ) - eJ - eX P^ 2(r 2 1 +a 2 ) )' 

A X (9) = e f -G+Hd6 = cos^m (nc^O^Holj 



where 



G 



H := V2rj A j3. (A10) 



2sin0(r| + a 2 )' 
The other holomorphic condition is 

= 3A a - fiXa = (6 + (3)\i - (i\ . (All) 

The solution has the form Ax = B_ x (8)e- i<l> / 2 + B 1 (8)e ict >/ 2 . 
It satisfies the differential equations 

<9 9 B_i - (G - H)B-i - V2fJ,fj A C-iA-i = 0, (A12) 
d e B 1 + (G + - V2/j,rj A C x A x = 0, (A13) 



where 



Bi{9) 



Ci 



V2A^ 



h(s)ds, 



(A14) 
(A15) 



and 



y h(9)d9 = J 2fj,r) A A 1 A_ 1 d9. (A16) 



USe ^ - ~^^+^ 2(ri+a 2 ) 2 ' We 8 et 



JJ'Ml^-i I 2 sin 6»d6» 



2ar(r 2 +a 2 ) 2 



(A17) 



and 



/ /i(s)<£s = / 2/xr? A ^i^-i^ 
Jo Jo 

= ^[ A /r^+ a 2: Q (2M-2r + l)-3r 3 ln( a + ^ 2 + a2 ) 



+ Vr 2 + ■ 



+2(r + M)a 2 ln(-a + V?' 2 + a 2 ) 



-Mr 2 In In ( 



a + \/r 2 + 



-a + V^ 2 + 



= ) - 2a 2 Mlna 



-2ra 2 ln(a + V?' 2 + a 2 )]/(a(r 2 + a 2 ) 2 ). 



(A18) 
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Therefore, we use —(A17\/ A18 to calculate Dougan-Mason 
mass. Thus we get 



and 



ttidm — 2Mr. 



/* -H(0) | A_i (0) | 2 sin 6d8 
\SZK0)M\ 



M- -^ + 0(a 4 ). 
2AM 1 ; 



J h(6)d9 = J 2nr) A A 1 A_- L d6, (A28) 

For Kerr solution in Bondi coordinate, we have \i — — ^A. 
Thus we get 



(A19) 



rriDM = 2Mr- 



3. Using the Dougan-Mason conditions to calculate the Kerr 
mass in the approximate Bondi frame 



The complex spatial null tetrad for Kerr in Bondi coordinate 
has the form 



M 



ft -Ai((9)|A_i((9)| 2 sin 6d6 
0( a 4 ) 



2AM 



(A29) 



V2r/A sm0 



which is approximate to the Kerr mass when a is very small. 



4. Checking our quasi-local formula 



(A20) 



The NP /3 is function of on horizon (recall eq. ( A4 1) and it 
has the form 

ofo . V2._, . \/2V A cot0 V2mcos(26») 

= -— iDsm9r] A + 



in order to check our 



We replace /i by using , 
quasi-local formula. Since we know the NP coefficients and 



Weyl tensors for Kerr in Bondi coordinate (See Section A 1 
we have 



4S A 



4sm0£ A 



sin 2 0E A D 2 



E 



Dougan-Mason's holomorphic conditions in components has 
two equations. One is 



24t/S A : 



= SA = (S - f3)\ . 



(A21) 



2e 3S A 24U 

where 

E := 8£ A - 4r A + 4a 2 (a 2 cos 2 - r A sin 2 0) . (A30) 



The solution has the form A = G_iA_i(0) e -^/ 2 ^ 
CiA l {6)e l ' t >/ 2 where A x ,A-i satisfy 

- (G + = 0, + (G - = 0. 

We use Maple to check the solutions are 

e jG+HdB 

sin(6»/2)^ A "exp(icos0D(3r A + a 2 cos 2 0)/12), 
cos (0/2) exp (i cos 0L> (3r A + a 2 cos 2 0) / 1 2) , 



A 1 {6) 



We thus get 



m,DM = 2Mr + 



M - 



j-tt _ t gl+a^+vrggg ] | (0) | 2 sin 0rf0 



fih(6)d0\ 



24M 



+ 0(a 4 ), 



(A31) 



G 



1 



2 sin 



(A22) 



The other is 

= 9Ai-mAo=(<* + /3)Ai-mV (A23) 

The solution has the form A x = J B_ 1 (0) e -^/ 2 + B^e^ 2 . 
Sx, B-i satisfy 

dgB-i - (G - H)B_i - s/2 mA C-xA- X = 0, (A24) 
+ (G + ff)Bi - \/2 Wa CiAi = 0, (A25) 



which is approximate to the Kerr mass when angular momen- 
tum per mass is very small. 

Remark. We conclude that the quasi-local mass for these 
three cases are the same up to a 2 . This is not surprising be- 
cause they are based on Dougan-Mason's quasi-local mass 
definition. Even though our work for the Kerr in Bondi coor- 
dinate is just approximation, we still can get the same quasi- 
local mass up to the second order. However, if one compare 
with the Hawking mass, the second order is different from the 
Dougan-Mason mass. Hawking mass of Kerr horizon is 



= M[ 



M - 



l + (l-a 2 /M 2 )s i 



8A/ 



+ 0(a 4 ). 



(A32) 



where 



BM = vSrJ hMds - 



(A26) 
(A27) 



APPENDIX B: QUASI-LOCAL MASS AND ENERGY FLUX 
FOR THE VAIDYA SOLUTION 



A good example to study the dynamical horizon is Vaidya 
solution. For this, we know that the dynamical horizon is 
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a space-like hypersurface if the black hole is growing. On 
the contrary, if black hole is contracting then it's a time-like 
hypersurface. In the later case the dominate energy con- 
dition will not hold so we restrict attention to the former 
case. The outgoing null normal of the Vaidya solution is 
£ a = 4- + (nj-n — v — 7"a) ^tt, we also have p=a=0 on 

horizon, and $nn = We now show that if the dominate 



<2(r'+r A ) 

and $ 00 = ^ 
energy condition is satisfied then f A 



> 0. From the Ray- 



chaudhuri equation C(Q(i) = —®fg\ — &abV + ^>ab^> — 
2e0(£) + ^R a bt L l h . The expansion of the outgoing null nor- 
mal of Vaidya is 0(« = —Rep 



'W 1 2 J 

r-2M(ii) ± r' 

vanishes on horizon. However, it's time derivative does not 



so it 



vanish. Here we have CeQ 



2rl ~ ' 



00 



< 0. If 



the dominate energy condition satisfied, then $oo > 0, i.e., 
ta > 0. Let V denote a vector field which is tangential to 
the dynamical horizon H + , is everywhere orthogonal to the 
foliation by marginally trapped surfaces and preserves this fo- 
liation. We can always choose the normalization of l a and 
n a such that i a n a — 1 and V a — £ a — fn a for some /. 
Since V a V a = —2/, it follows that H + is space-like, null 
or time-like, depending on whether / is positive, zero or neg- 
ative. We now show that if the dominate energy condition 
holds then / is non-negative. Let us begin by noting that the 
definition of V a immediately implies Cv®(t) — 0, whence, 
C(Q(tj = /Ai@M' Therefore, the Raychaudhuri equation 
for £ a implies fC n ®(t) = H CT I 2 — $00 < 0. Physically, we 
can assume that the expansion of the outgoing null normal be- 
comes negative if it moves along the incoming direction to the 
interior of the marginally trapping surfaces, so C n Q^ < 0. 
Hence, if the dominate energy condition is satisfied, then / is 
non-negative, i.e., the dynamical horizon H + is space-like or 
null. The Vaidya solution is similar with the Schwarzschild 
solution but the mass is time dependent and it's a non-vacuum 
dust solution. In this appendix we calculate our quasi-local 
mass and flux for the Vaidya solution and compare with the re- 
sults one expects from the form of Vaidya metric. In advanced 
Eddington-Finkelstein coordinates (v, r, 0, <f), the metric is 

ds 2 = _ 2M(vl) dv 2 _ 2dvdr - r 2 (d6 2 + sin 2 9d(j) 2 ). 

(Bl) 

where the horizon radius r A (v) — 2M(v) is time dependent. 
When r approaches the horizon radius r A (v), the metric ap- 
proaches to 



ds 2 = - 2r A dv 2 - r\(d6 2 + sin 2 9d(j) 2 ) 



(B2) 



where we use a coordinate transformation dr' = dr — r A dv. 
It's obvious that the Vaidya horizon is space-like or null hy- 
persurface since it has the metric signature ( ) if r A > 0. 

It's null when the mass approach to constant and it returns to 
the Schwarzschild solution. In the new coordinate (v, r' , 8, </>), 
the tetrad components are 



t a = 



n 

m° 



V x ' 2(r'+r A ) 

(o , -1 

(0 , 



r A 






1 



0) 

0) 



The NP coefficients are 



\/2(r'+r A ) sin 

(B3) 



_ r-2M(v) _ 
P ~ 2r 2 — 1r 



/' 

e = 



1 . 

r 

M(v) 
2r 



1 

4r A 



O(r'), 
e - + 0(r'), 

^ + 0(r'). 



7 = 0, 
+ 0(r' 3 ), 

I2\ 



*0 = 
$00 



(B4) 



M(v) 

1 A 



O(r'), 



Quasi-local mass of Vaidya horizon 

Let's calculate the horizon energy momentum near the Vaidya 
horizon and examine the related mass gain or loss formula 
near such dynamical horizon. We use the holomorphic gauge 
fixing, so the horizon quasi-local energy-momentum is 



p AA' _ _ 1 



/iX^X^rdS. 



4tt j s 

Here we use the Dougan-Mason definition of mass: 



(B5) 



m 2 DM = P AA ,P^ = 



\N\ 



_( P Q3Lp 



11' 



pSLpiS-) (B6) 



where N = e— = XqXj — XjX^. To calculate the Dougan- 
Mason mass, we have to solve the holomorphic equations for 
the constant spinors equations. Here we simply follow the 
calculation that done by Bergvist J3). The holomorphic con- 
ditions: 



,0,1 



5A + 0-A1 = 
5Ai — (iXq = 



(S-/3)X = 
(5 + /3)Ai-^Ao = 



For 5Ao = 0, we have to solve the PDE 



d 



36 sin^, 







cot ( 



)A Q = 0, 



(B7) 
(B8) 



(B9) 



then we get the solution Ao = C_iA[j + C\X\. For 5Ai 
jLtAo = 0, we have to solve the PDE 



-[( 







d 



cot ( 



)X 1 + V2X Q ] = 0, (B10) 



y/2r A 86 sin 6 d<j) 2 
then we get Ai = C-iXj + C\X\. Here 
Xl = A^e-'^ 2 , xl = A x e^ 

,0 _ J e B(s)ds e-^ /2 ■ 1 _ fl B(s)ds 

A l ~ V2A 1 (6) ' A i ~ V2A^(6) ' (Bll) 

At = cos(6»/2)/r, A-i = sin(6>/2)/r, 

B{6) = sin 9/r 2 . 

We choose two independent solution C_i = 1, C\ — for 
A~ and C_i = 0, C\ = 1 for A~ here. So when we integrate 
over (/>, <j) < (j> < 2tv in P^-, we find = P^- and 
pQ0_ _ pJJ_ Therefore, the Dougan-Mason mass becomes 



i) 



m DM 



V2P 



no 



\N\ 



(B12) 
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where 



pOT/ _ _ 



■h I ^Wik (B13) 



and N = A§Ay - Af A5 = ^ Jfi h(9)d6 = -^J. Finally, the 
Vaidya horizon mass from the Dougan-Mason's mass defini- 
tion is 



m DM 



= M(v) 



(B14) 



which agrees with the expected result. 
Quasi-local flux of Vaidya horizon 

Next, we reduce our flux formula to Vaidya by taking shear 
vanishing and compare it by using the solution of the holo- 
morphic constant spinors. For Vaidya solution, we have 



1 



Mo 



r A (v) 



(B15) 



where r A is time dependent. We substitute these formulae for 
the NP quantities from the Vaidya solution into our quasi-local 
energy-momentum integral we obtain in (|68j). Then 



1 

47T 



A i oAoA ,(iS'A 



(B16) 
(B17) 



Then we rescale A{j so that we get the condition Aq = (recall 
d54]l) from our asymptotic constant spinor condition for the 



dynamical horizon. We use the fact that j^dS A — 2^-dS A . 



From the first expression (B16 1, the quasi-local energy flux is 



dl 
dv 



~ /(Ao + 2mo — )AgA°,d5 A (B18) 
47r J r A 



(B19) 



From the second expression (B 17 1, it becomes 



dl 

dv 



1 

47T 



2e 2e 2 Q 



^ ■ ^)A%dS A 



(Use = 
1 /" M 



8M 2 



eo^A 
1 

8M" 1 



4tt / 2M 2 



AgAQ/dS'A- 



We use the Dougan-Mason's mass definition, then tudm = 
^pfc . Therefore we can put A§ = A§ in our quasi-local 
formula. The Vaidya's quasi-local mass flux is 



dm 



DM 



dv 



r 

47T 



M sm(d/2f 
2M2 r 2 



dS = M. (B20) 



If the dominate energy condition satisfied, i.e., M > 0., then 
we find the mass gain law for the Vaidya. 
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